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Abstract 



The problem of assigning a momentum to an electromagnetic wave 
packet propagating inside an insulator has become known under the 
name of the Abraham-Minkowski controversy. In the present paper we 
re-examine the question, first through a power expansion in the polar- 
izability of the medium and assuming the simplest and most natural 
choice for the force exerted on a dielectric material by an electromag- 
<— j | netic field. It is shown that the Abraham expression is highly favoured. 

r Q-f We then show the complete generality of these results. 

1 Introduction 

It is well known that an electromagnetic (e.m.) wave propagating in 
the vacuum carries both energy and momentum jl], [2]. Of course this 
must be true also for radiation propagating inside matter. Soon af- 
ter the formulation of Einstein Relativity Theory, conflicting expres- 
sions were proposed for the space momentum to be assigned in this 
situation [3], jl]. This discrepancy, known as the Abraham-Minkowski 
controversy, is still attracting a considerable interest in part of the 
scientific community^ 

In order to discuss the problem we imagine an idealized situation in 
which matter is schematized as a uniform, isotropic dielectric material 



e-mail: massimo.testaQromal .infn.it 
2 For a review of the subject see[5],[6] and [7]. 
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T>, occupying the space region Vp, with dielectric constant e, refraction 
index n = y/e, with no magnetic properties, i.e. with permeability 
/j, = 1 and sufficiently massive that it can transfer space momentum, 
but not energy with the incident radiation. The dielectric is surrounded 
by empty space, V v . If stress variations induced by e.m. radiation 
can be safely neglected, the Abraham- Minkowski discrepancy can be 
formulated imagining a collimated e.m. wave packet, coming from the 
vacuum, V, with an energy u and a momentum p (u = c|p|), entering 
the dielectric material. Once this radiation is fully inside the dielectric 

• Minkowski attributes to it the same vacuum energy u and an 
absolute value of the momentum n|p|, 

• while Abraham attributes to it again the same u, but and an 
absolute value of the momentum — . 

n 

It is the purpose of the present paper to give a small contribution 
to this debate through the discussion of examples, amenable to quanti- 
tative conclusions, which clearly favour the Abraham point of view on 
the momentum. Approximate computations are performed in sections 
13. 1^ 13.21 13. 3[ under the assumption that the dielectric polarizability of 
the medium is small. In section H] the problem is discussed without 
approximations, studying the transition of the e.m. field through the 
dielectric boundary in full generality. 

2 Polarization charges and currents 

The description of a transparent, dielectric, homogeneous material with 
dielectric constant e and magnetic permeability fi = 1 requires the 
introduction of the polarization field, P(x, t), in terms of which the 
polarization charge density pp(x, t) is given as 

p P (x,t) = -V-P(x,t). (1) 
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The time dependence of the polarization field gives rise, as a conse- 
quence of charge conservation, to a volume current density described 
by 

J P (x,t) = P(x,t). (2) 
The Maxwell equations in presence of polarization are written as 

V ■ E = -V • P V ■ B = (3) 

VxE = -±B VxB = -E + -P. (4) 

c c c 

In the linear regime, P(x, t) is related to the electric field E(x, t) 
through 

P(x,t) = (e-l)E(x,t) =a E(x,t), (5) 

which defines the polarizability a. 

In the following we will consider a uniform dielectric, i.e. we will 
consider the polarizability a as constant, except in a thin region around 
the dielectric boundary, in which it quickly drops from a to 0. When- 
ever possible, without giving rise to any ambiguity, we will treat this 
transition region according to distribution theory [T3] and will schema- 
tize the situation through the introduction of a surface density op of 
polarization charges induced on the dielectric surface E 

p P (x, t) = -V ■ P(x, t) = a (n v ■ E(x, t)) «J(E(r)) = cr P (x, t) 5(E(r)) , 

(6) 

where is the outward normal to the dielectric surface, E(r) = 0, 
and 5(E(r)) is the surface Dirac 5-function[13j with support on S. In 
eq.(j6]) we defined 

ap(x,t) =a K-E(x,t)). (7) 

The electric field E(x, t) appearing in eqs.fjS]) and (J7|) is thought to be 
obtained reaching the surface E from the dielectric side. 
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According to eq.(j6]), polarization charges are not present in the 
dielectric interior 

p P (x, t) = -V • P(x, *) = -a V • E(x, t) = 0, (8) 

in virtue of the hypothesis that the sources generating the incident 
electromagnetic field are external to the dielectric [9]. 

We remind the reader that across the dielectric boundary we have 
the following jump conditions [TJ 

AE n = a P AE|| = (9) 
AB = 0, (10) 

where E n and E|| denote the normal and tangential components of the 
electric field on E and the symbol A denotes the increment in passing 
from T> to V. 

3 The forces on a dielectric material and the a- 
expansion 

One fundamental issue raised in the literature, concerns the precise 
form of the e.m. force acting on a material, in terms of the matter 
polarization [8] , [9] , , [EE] , [Lj . 

We will take the simplest possible attitude, assuming that the e.m. 
field acts on the polarization charges and currents in the same way 
as it acts on the so called free charges and currents[14j. Moreover, 
as discussed in the introduction, we neglect forces due to a change of 
the stress status of the dielectric body. The total electromagnetic force 
acting on the dielectric, occupying the region Vp, is therefore computed 
as 

F v (t) = f p P (x,t)E(x,t)dr+- / Jp(x,t) x B(x,t)dr = (11) 

JVd C JV-d 
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[E • Va] E dr + 




ExBdr. (12) 



If we consider the transition from the dielectric to the vacuum as sharp, 
i.e. if we insist using eq.flH]), we must consider eqs.(!TTT) and (T12]) as for- 
mal. In fact pp(x, t) is a distribution singular on E, as shown in eq.(j6]) 
or, equivalently, the term Va appearing in eq.f lT2"j) . is a distribution 
singular on S. The problem is that in eqs.f llip and ( I12p these distri- 
butions appear multiplied by E, which is itself a distribution singula^ 
on S. This is an illegal operation in distribution theory and could give 
rise to ambiguities or infinities. In the present case, for example, it is 
not clear if the E(x, t) appearing in eqs.ffTTT) and ( !T2|) should be taken 
from the vacuum or the insulator side. 

We will not need, at the moment, the full solution to this problem, 
which we postpone to section HI In fact, in order to be able to perform 
computations, we will introduce an approximation scheme. The tool 
we will use is to imagine an homogeneous, weak dielectric, namely a 
material with a polarizability such that 



Eq.(fl3l) allows us to consider an expansion in powers of a, first used by 
J. P. Gordon[15j. We will apply this approximation scheme up to first 
order in a, which is easy and accurate enough to distinguish between 
the Abraham or Minkowski formulation and is not plagued by any 
ambiguity, as it will be clear in a moment. 

We consider therefore a dielectric block invested by an e.m. wave 
described by a given electric field E (x, t) and a given magnetic field 
B (x, t), generated by sources external to the dielectric itself. 

Since we deal with infinitely massive matter, it is not possible to 
trace directly the momentum flow coming from radiation. We can, 

E is discontinuous across E, see eq.©. 



a = e — 1 



0. 



(13) 
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instead, identify the transferred momentum through the knowledge of 
the force that the e.m. wave exerts on the dielectric |1 6]. |1 7]. 

Up to order a, the force acting on the dielectric T>, formally given 
by eq. f|T2|) . reads 

F v (t)m-f [E ■ Va] E dv + - [ E xB dr = 

JVt, C JVt, 

= a f(E -n v ) E cE + - / E xB dr, (14) 

•/£ C JVt, 

which shows that, up to order a, no ambiguity is present, since E is 
perfectly regular on E. 

The general treatment of the singularity problem, present in eqs.f llip 
and (11 2p will be given in section HI 

The use the identity 

- E x B = (E • V)E - iv(Eg) , (15) 

which follows from the source-free Maxwell equations, allows us to com- 
plete the time derivative in eq. ([T4"]) and write Fjj(t) as 

E xB rfr + ^/ V(E^)dr= (16) 
c at JVt, 2 Jvt> 

= --/ E xB dr+- n P E (IS. 17 

c dt Jv-d 2 Jt, 

In the following sections we will use eqs. (EES)) and (Fl7]l to analyze simple 
situations, obtaining results in consistent agreement with the Abraham 
form of the momentum. 

3.1 The Einstein box 

In a beautiful paper, Balasz[T8] considered a thought experiment which 
consists in sending an e.m. wave packet on a dielectric box. The input 
used by Balasz in order to predict the outcome, is the requirement that 
the combined center of mass of the box and of the radiation should 
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maintain, during the experiment, a uniform, rectilinear motion. On 
this basis Balasz concluded that, while the radiation penetrates into 
its interior, the box should be subject to a force in the direction of the 
e.m. wave, it should then not be subject to any force while the wave 
is traveling inside the dielectric and should finally be subject to an 
equal and opposite force, during the time in which the radiation exits 
from the box. We will analyze this problem from the point of view 
of eq.f ll6p . finding results in agreement with Balasz analysis and the 
Abraham form of the moment uuQ. 

We consider an unperturbed e.m. wave packet propagating along 
the positive z-axis, of the form@ 

E (x,*) = (/(2-ct),0,0) (18) 
B„(x,*) = (0,/(z-cO,0), (19) 

where the function / has a finite support 

f(z-ct)^0, a<z-ct<b. (20) 

Imagining Vp as a cube of side L, with two faces orthogonal to the 
z-axis, we have, from eq.( fT6l) 

FvM~-tJ 9(z-ct)dV+- —g(z-ct)dV = 
c dt Jv-v 2 Jv-p oz 

a 



where g(z — ct) = f (z — ct), so that 



/ g'(z-ct)dV, 



F Vz (t) w a ^ [g(-ct) - g(L - ct)] , (21) 

where A, smaller than L 2 , is the transverse section of the wave packet. 
If in eq. (l2Tjl we take L > b — a, then g(—ct) and g(L — ct) cannot 



4 In this section we follow ref.|15| 
In the spirit of the a-expansion, Eo(x, t) and Bo(x, i) are the electric and magnetic 
fields at zeroth order, a , and therefore they freely propagate at the speed of light in 
empty space, for all t. 
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be both different from for any t and we can consider three distinct 
regions of time: 

• Region I 

— - < t < — -, during which the e.m. wave is entering in Vp and 
a positive force is exerted on the dielectric, 

• Region II 

— - < t < — -, during which the e.m. wave is traveling inside Vv 
and no force is exerted on the dielectric and 

• Region III 

< t < during which the e.m. wave is exiting from Vp 
and a negative force is exerted on the dielectric, 

in agreement with the Balasz analysis. 

3.2 The e.m. momentum inside the Einstein box 

In the Einstein box experiment discussed in section l3TTj the momentum 
of the incident e.m. wave packet is directed along the z-axis and its z- 
component is given by 

= - f g{z - ct) dV = - f g{z)dz . (22) 

' C JV-d C J a 

During the time region I, the momentum transferred to the dielectric 
block through e.m. forces is also directed along the z direction and is 
given by 

p v = n Fv ^ t)dt =a \ n a{ - ct)dt = i p ^ o) • (23) 

c c 

Therefore, if we require momentum conservation, we are led to attribute 
to the e.m. wave, while it propagates inside the dielectric during the 
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time region II and not exerting any force on it, a ^-component of the 
momentum such that 

P (ii) = p(o) _ p <p = p(o )(1 - |) « ^ . (24) 



In eq. p4jl we used the relation 

n = yfl « 1 + - . (25) 

relating the refraction index n to the dielectric constant and we recov- 
ered the Abraham form of the radiation momentum. 

3.3 The case of oblique incidence 

In this section we discuss, on the grounds of a more general kinematic 
configuration, the validity, up to order a, of the Abraham form of the 
radiation momentum. For this purpose we consider an e.m. wave packet 
incident on a dielectric boundary in an oblique direction. The dielectric 
occupies the half space z > 0, always denoted by Vp, and, being as 
before infinitely massive, it does not exchange energy with the e.m 
wave hitting ito- The interest of this more general kinematical situation 
resides in the fact that it allows to recover, besides the Abraham form 
of the momentum, also the Snell's law for refracted light. 

We start integrating eq.flTTj) over time, from — oo to +oc0. We get 



/+oo 
F v (t)dt 
-oo 



a 



E x B \ t=+oo dr- E X BqI^ dr) + 
+n v - dt / Eg dS = 
- E x B |, =+oo dr-z- dt / E^dS (26) 

C JVn ^ J —oo JYi 



6 This point will be further discussed in section [4] 

7 The t-integration is actually limited to the time interval during which the wave packet 
crosses the boundary z = 0. 
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where Px> is the momentum transferred to the dielectric by the radi- 
ation, S is the dielectric surface, z — 0, np its external normal and 
z = — ri£>, the unit vector along the positive 2-axis. Since V© is the 
region occupied by the dielectric, the term J v E x Bo| t= _ OG dr is zero, 
because, for sufficiently negative t, the radiation is supposed to be con- 
fined in the region z < 0. Similarly we have 

~ f E xB | i=+oo rfr = p(°) ; (27) 

where P^ -* is the initial e.m. momentum. In fact as the time evolves 
we imagine that the unperturbed e.m. radiation will, eventually, be en- 
tirely inside Vr>: since the electric and magnetic fields E (x, t), B (x, t) 
propagate freely, as explained in the note El the integral appearing in 
eq. (1271) coincides with the initial, vacuum momentum of the unper- 
turbed e.m. field. 

Summing up we have 

Pp = «Pf -z-/ dt EjjdE. (28) 

2 J —oo JTi 

In order to give an estimate of the last term in eq. (!28|) . we use the 
approximation that the wave packet propagates without deformation, 
while crossing the dielectric boundary. In other words we assume 

Ejj(x,t) = D(x- ctk) , (29) 

where k is the unit vector which defines the (incident) e.m. wave di- 
rection. 

Eq. (I29p is valid for wave packets containing short wave lenghts and 
propagating over a not too long time-lapse and it will be a good ap- 
proximation for e.m. wave packets with small space extension. 

Performing the change of integration variables y = x — ctk, we 
easily get 

/+oo p r+oo r 

dt E 2 dZ= dt D(x - ctk) dZ = 
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,+oo 77(0) p(0) 

, , dyD(y) = ^- = ^-, (30) 

where L^ -* = cP^ is the e.m. energy of the free initial wave packet 
and 

k z = cos i , (31) 

i being the angle between the direction of the incident e.m. wave and 
the z-axis. Therefore we obtain, for the total momentum transferred 
to the dielectric wall, 

P^aPf-zApf. (32) 
2 cos i 

By the requirement of momentum conservation, we are therefore led to 
attribute to the e.m. wave, while inside the dielectric, a momentum 
such that 

p(0) = P' + Pv = P ; + a P (o) _ z -^p(°) , (33) 

2cosi 

so that 

Pi = (1 - a) Pf + z -^Pf . (34) 

2 cos i 

Let us consider the implications eq. (1341) . First of all we have, to order 



(P;) 2 = (1 - 2a) (Pf) 2 + ^ P W(P(% = 

COS I 

1 - 2a) (P<°>) 2 + a(Pf) 2 = (1 - a) (Pf ) 2 = -J^f , (35) 



n- 

which again confirms the Abraham expression for the radiation mo- 
mentum. Moreover, from eq. (l34|) . we also get 

,, (K)± (l-«)(PfU 



sin«' - 



PI, " (l-f)(Pf) 

a, (Pf)i l(PfU 1 . , ,„ R , 

1 ) — = — = — smi, (3o) 

v 2 ; (pf) n (pf) n 1 ; 
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where i' is the refraction angle. Eq. fl36|) shows the validity of Snell's 
law, starting from the conservation of the Abraham momentunj^. 



4 Exact results 

Previous sections suggest that, up to order a, once the e.m. wave crosses 
the dielectric boundary, its associated momentum changes from the ini- 
tial momentum p the wave possessed in the vacuum, to the Abraham 
momentum value inside the dielectric. This has been checked indi- 
rectly, assuming momentum conservation and attributing to the e.m. 
wave the missing momentum. The physical picture associated with 
the change in momentum of the e.m. wave, is that the incident e.m. 
radiation, once it reaches the dielectric boundary, excites polarization 
charges and currents which induce an e.m. field carrying the additional 
missing momentum. 

In order to study this problem we should approximate the e.m field 
up to order a 



E ~ Eq 4~ Ei 



(37) 



BwB + B 




This can be done [19] and also the Abraham Poynting vector 



n = -E x B, 



(39) 



c 



should be expanded according to 



n « n + n 



(40) 



where 



o — 



- E n x B 



(41) 



c 



In section[l]we show that, up to order a, there is no reflection. 
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and 

n x = -Ej x B + -E x Bi. (42) 

c c 

Eq.( l42l shows that, up to order a, there cannot be any reflected wave: 
the support of IIi is contained in that of n and therefore the support 
of the momentum density II, coincides with the support of the zeroth 
order e.m. field, Eo and Bo at all times. 

Although the small a-approximation is very useful to get a semi 
quantitative idea of the basic phenomena related to dielectric polariza- 
tion, it is not easy to establish the degree of generality of the results 
obtained in this way. 

In this section we will overcome the limitations due to the a-expansion 
and discuss in general terms the passage of an e.m. wave packet from 
vacuum to an insulator. 

Following Abraham^ we make the hypothesis that the e.m. wave 
momentum is given in any situatio by J noc i Ex B dr, where 
denotes full three dimensional space and we will trace the momentum 
change while e.m. wave crosses the boundary between the empty region 
Vy and that occupied by the dielectric, Vp- In this way we will be able 
to give a meaning to the formal expression for the force exerted on 
T> given in eqs. fllll) and (|12|) . obtaining an unambiguous result, which 
coincides with that found in texbooks|20| in similar situations. 

In order to deal with problems related to the discontinuity of the 
electric field across S, we introduce a regularization. In other words we 
isolate the surface S enclosing it inside a thin layer £$ with thickness 
5 and will eventually take the limit 5 — > 0, in which £$ tends to E. 

We denote by IZoo Cs the full three dimensional space TZoo, de- 
prived of C$ and consider the quantity J^q/v \ E x B dr. If the inte- 
gration region were TZoo, it would represent the Abraham momentum 



9 It is easy to use the same procedure with the Minkowski momentum. 
10 We are considering configurations in which the e.m. field is confined in a finite region. 
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of the e.m. wave. It is clear, at least from physical intuition, that 

lim / - E x B dr = ; (43) 
s-yoJc s c 

in other words we do not expect a finite fraction of the e.m. momentum 
to be trapped and accumulate inside £5 with a vanishing 5. This means 
that lim^o S-jiooQCs E x B <ir correctly gives the Abraham momentum 
of the e.m. field. 

We are now in the position to compute 

d r 1 „ , 

— / - E x B dr = 

dt Jiz^eCs c 

= [ -ExBdri / -ExBdr, (44) 

JllovGCs c J-RcoeCs c 

where E and B are perfectly regular, having isolated the regions con- 
taining their discontinuities. 

Through the use of the Maxwell eqs. and (jlj) we transform some 
of the volume integrals into surface integrals and get 

d r 1 
lim — / - E x B dr = (45) 
s^o dt Jn^eCg c 

= f l-l-n v E 2 v + [(n v ■ E v ) E v ] dH + {V V} + 

ExBdr, 

where Ex> and are the electric field and the external normal respec- 
tively, computed on £ as a limit from the dielectric side and {T> — > V} 
represents the same expression, computed as a limit from the vacuum 
side. There is no term containing the magnetic field, as a consequence 
of its continuity across S, eq. pUj) . 
Using eqs. §3§ we get 

d r 1 
lim — / - E xBdr = 

s^o dt Jn^eCg c 
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/[-in p AE 2 + A(£ n E)]d£- - — -/ ExBdr = 

is 2 c 

-- f a P (E v + E v )dZ-- [ 3 P xBdr = -F v (t), (46) 



where 

l r _ , „ l 



F©(t) = 77 / a P (E v + E^dE + - / JpxBdr (47) 

2 JS C JVb 

= - / dp (Ey + E c )rfS + / ExBdr. (48) 

2 je c 7Vu 



We now imagine that at an initial time T\ the radiation is entirely con- 
tained inside the vacuum region Vy and at a final time T 2 the radiation 
is entirely contained inside the dielectric region Vx>. We are in this 
way excluding any substantial radiation reflection from the dielectric 
boundary. While, as we have shown, up to order a this condition is 
automatically satisfied, this is not assured for the exact solution and 
we must explicitly assume that the experiment is performed under neg- 
ligible reflectivity conditions. 

Integration of eq.f )46p over t from T\ to T 2 gives 

lim / -Ex B|r 2 dr — lim / -Ex B| T dr = 

5->o JKaceCs c s^oJn^eCg c 

rT 2 



/ 2 F v (t)dt (49) 



which is interpreted as 



P 7 ' - P 7 = - r F v (t)dt = -P v , (50) 

suggesting the identification of Pp with the momentum transferred to 
T>. This result makes eq. (1471) a very natural candidate for the force 
exerted by the e.m. field on an insulator. Eq. (l4"71) provides the re- 
quired regularization of the force acting on the dielectric surface S: it 
must be computed using the average of the electric field on E from the 
vacuum and from the dielectric side. This is the natural choice when 
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computing the force acting on a simple layer charge distribution, as 
discussed in textbooks [20j. With this prescription eq.flUJ) represents 
the precise statement that polarization charges and currents interact 
in an universal way with the e.m. field. 

As a consistency check we notice that Fx> = as soon as the e.m. 
field is completely internal to Vp. This is obvious for the first term of 
eq. fHHj) . As for the second term, J v E x B dr, when E and B do not 
touch the boundarj0 E, can be transformed, using Maxwell equations, 
as 

/ ExBdr=- [ (VxB)xBdr = 0, (51) 

Jv-v e Jv-r, 

through partial integration. Therefore, as soon as the e.m. wave aban- 
dons S, it ceases to exert any force on the dielectric body, in agreement 
with what was found in the approximate computation of section 13.11 

These arguments are easily extended to the study of the correspond- 
ing energy flow. 

The e.m. energy density, U, for a dielectric material is given by[2] 

W = ^(eE 2 + B 2 ). (52) 
In analogy with eq.flUJ), we compute 

4- - [ (eE 2 + B 2 ) dr = [ (eE ■ E + B ■ B dr = 
= - f [(E v x Bp) • n v + (Ey x B v ) • n v ] = 

= - f ME x B) • nj>] dT, . (53) 
As a consequence of eqs.flH]) and ffTUl) we have 

A[(E x B) • tl v ] = (AE x B) • ri£> = a(n v x B) ■ n v = , (54) 



11 The important point is that e must not show discontinuities within the region of 
integration. 
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so that 

limj / (eE 2 + B 2 )| T2 rfr = limJ / (eE 2 + B 2 )| Tl dv , (55) 

which shows that, again under the assumption of negligible reflection, 
there are no energy losses in the radiation during the transition from 
V to V. 

5 Conclusions 

We presented some elementary checks supporting the validity of the 
Abraham expression of the momentum to be attributed to an e.m. 
wave packet, in order to obey momentum conservation. We used the 
approximation e ~ 1, first proposed in ref.[15], which allowed the ana- 
lytic treatment of some problems. In this way it has been possible to 
derive the Snell's refraction law through momentum conservation. 

In order to overcome the limitations due to the expansion in powers 
of the polarizability, in section H] we presented a general, nonperturba- 
tive argument showing that Abraham momentum conservation implies 
that an e.m. field interacts with polarization charges and currents in 
the same way as it does with the so called free ones. The argument 
shows the correct way of dealing with the discontinuity of the electric 
field across the dielectric boundary. 
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